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SECTION  I 


INTRODUCTION 

The  original  Rotor-Bearing  Dynamics  Design  Technology  Series  AFAPL-TR-65-45 
(Parts  I  through  X)  included  a  volume,  Part  IV(l)f  which  presented  design 
data  for  typical  deep-groove  and  angular  contact  ball  bearings.  The  data 
was  presented  in  graphical  form  and  consisted  of  direct  radial  stiffness, 
load  carrying  capacity,  and  load  levels.  In  addition  design  guidelines  and 
limitations  were  discussed.  The  major  deficiencies  of  this  original 
volume  were  that  centrifugal  effects  due  to  high  speed  were  Ignored,  and 
axial  and  angular  stiffness  Information  were  omitted. 

Subsequent  to  the  publication  of  Part  IV,  several  extensive  treatments  of 
rolling  element  bearings  including  elastohydrodynamic,  thermal,  and  cage 
effects  have  been  published.  The  computer  program  of  Mauriello,  LaGasse, 
and  Jones  (2)  considers  both  elastohydrodynamic  and  cage  effects  for  ball 
bearings.  The  more  recent  computer  based  design  guide  prepared  by  Crecelius 
and  Pirvics  (3)  treats  elastohydrodynamic,  thermal,  and  cage  effects  for  a 
system  of  ball  and  roller  bearings. 


(1)  Lewis,  P.  and  Malanoski,  S.B.,  "Rotor  Bearing  Dynamics  Design 
Technology.  Part  IV:  Ball  Bearing  Design  Data  Technical  Report," 
AFAPL-TR-65-45,  Part  IV.  Air  Force  Aero  Propulsion  Laboratory, 

Wright  Patterson  AFB,  Ohio. 

(2)  Mauriello,  J.A.,  LaGasse,  and  Jones,  A.B.,  "Rolling  Element  Bearing 
Retainer  Analysis,"  DAAJ02-69-C-0080 ,  TR105.7.10,  USAAMRDL-TR-72-45. 

(3)  Crecelius,  W.T.  and  Pirvics,  J.,  "Computer  Program  Operation  Manual 
on  "SHABERTH"  a  Computer  Program  for  the  Analysis  of  the  Steady  State 
and  Transient  Thermal  Performance  of  Shaft  Bearing  Systems,"  AFAPL- 
TR-76-90,  Air  Force  Aero  Propulsion  Laboratory,  Wright  Patterson  AFB, 
Ohio,  October  1976. 


Thus,  very  sophisticated  analytical  tools  are  available  for  the  design  and 
application  of  rolling  element  bearings.  Neither  of  these  tools,  however, 
provide  the  user  with  the  stiffness  matrix  required  for  solution  of  rotor 
dynamics  problems.  In  addition,  both  computer  programs  are  very  large  and 
require  an  extensive  computer  facility  for  use. 

Part  11(4)  of  the  revised  series  provided  an  update  of  the  original  Part  IV 
(1).  Those  aspects  of  the  original  Part  IV(1)  which  treated  general  design 
aspects  of  ball  bearings,  load  capacity,  speed  limitations,  etc.  were  deleted 
since  their  coverage  is  superficial  compared  to  the  more  sophisticated  com¬ 
puter  tools  now  available  (2,3).  Only  those  parts  directly  connected  with 
preparation  of  input  for  the  rotordynamic  response  programs  (Part  1(5)  of 
the  revised  series)  were  retained.  The  stiffness  data  included  in  the 
original  Part  IV  were  also  updated.  A  later  volume  (Part  111(6)  of  the 
revised  series),  enlarged  the  treatment  to  include  the  tapered  roller  bearing. 

The  present  volume  (Part  IV  of  the  revised  series)  extends  the  treatment  of 
rolling  element  bearings  to  the  cylindrical  roller  bearing.  The  complete 
stiffness  matrix  is  calculated  including  centrifugal  effects.  Considerations 


(4)  Jones,  A.B.,  and  McGrew,  J.M.,  "Rotor  Bearing  Dynamics  Technology  Design 
Guide — Part  II:  Ball  Bearings,  "AFAPL -TU-73- 6,  Part  II,  February 

Air  Force  Aero  Propulsion  Laboratory,  Wright  Patterson  Air  Force  Base, 
Ohio. 

(5)  Pan,  C.H.T.,  Wu,  E.R.,  and  Krauter,  A. I.,  "Rotor  Bearing  Dynamics 
Technology  Design  Guide:  Part  1,  Flexible  Rotor  Dynamics,"  AFAPL-TR- 
78-6,  Part  I,  June  1978,  Air  Force  Aero  Propulsion  Laboratory,  Wright 
Patterson  Air  Force  Base,  Ohio. 

(6)  Jones,  A.B.,  and  McGrew,  J.M.,  "Rotor  Bearing  Dynamics  Technology  Design 
Guide  -  Part  III:  Tapered  Roller  Bearings,"  AFAPL-TR-78-6 ,  Part  III, 
February  1979,  Air  Force  Aero  Propulsion  Laboratory,  Wright  Patterson 
Air  Force  Base,  Ohio. 


such  as  elastohydrodynamic  and  cage  effects  are  not  included  since  they  have 
little  influence  on  the  calculation  of  cylindrical  roller  bearing  stiffness. 
The  resulting  program  (Appendix)  is  reasonably  small  and  easy  to  use. 


SECTION  II 


ANALYSIS 

2.1  General  Bearing  Model  and  Coordinate  System 

Accurate  calculation  of  the  lateral  dynamic  response  of  a  high-speed  rotor 
depends  on  realistic  characterization  of  the  support  bearings.  In  the  most 
general  case,  both  linear  and  angular  motions  are  restrained  by  the  support 
bearings  at  the  attachment  location.  In  the  analytical  model,  the  reaction 
force  and  the  reaction  moment  of  each  bearing  are  felt  by  the  rotor  through 
a  single  station  of  the  rotor  axis.  As  schematically  illustrated  in  Figure 
la,  a  coil  spring  restraining  the  lateral  displacement  and  a  torsion  spring 
which  tends  to  oppose  an  inclination  are  attached  to  the  same  point  of  the 
rotor  axis.  A  complete  description  of  the  characteristics  of  the  support 
bearings,  however,  involves  much  more  than  the  specification  of  the  two 
spring  constants.  This  is  because: 

.  The  lateral  motion  of  the  rotor  axis  is  concerned  with  two 
displacement  components  and  two  inclination  components. 

.  The  restraining  characteristics  may  include  cross  coupling 
among  various  displacement/inclination  coordinates. 

.  The  restraining  force/moment  may  not  be  temporally  in  phase 
with  the  displacement/inclination. 

.  The  restraining  characteristics  of  the  bearing  may  be  dependent 
on  either  the  rotor  speed  or  the  frequency  of  vibration,  or  both. 

.  Bearing  pedestal  compliance  may  not  be  negligible. 

To  accommodate  the  above  considerations,  the  support  bearing  characteristics 
are  described  in  Reference  5  by  a  four-degrees-of-freedom  impedance  matrix 
as  defined  in  Equation  (1): 

’  14,  a) 


-4- 


Figure  la  Bearing  Stiffness  Model 


Figure  lb  Bearing  Location  Coordinate 


where  is  a  column  vector  containln3  elements  which  are  the  two  lateral 

displacements  (6,6)  and  the  two  lateral  inclinations  (0  ,  6  )  of  the 
x  y  x  y 

rotor  axis  at  the  bearing  station  N. 


Employing  a  right-handed  Cartesian  representation  in  a  lateral  plane  as 

depicted  in  Figure  lb,  the  z-axis  is  coincident  with  the  spin  vector  of  the 

rotor.  The  x-axis  is  oriented  in  the  direction  of  the  external  static  load, 

and  the  y-axis  is  perpendicular  to  both  z  and  x  axes  forming  the  right-handed 

triad  (x,  y,  z) .  (6,6)  are  respectively  lateral  lineal  displacement 

x  y 

components  of  the  rotor  axis  along  the(x,  y)  directions.  (0,0)  are 

x  y 

lateral  inclination  components  respectively  in  the  (z-x,  z-y)  planes.  Note 

that  0  is  a  rotation  about  the  y-axis,  while  0  is  a  rotation  about  the 

y  }  x 

negative  x-axis. 


ZN  is  a  complex  (4x4  matrix) ,  and  in  accordance  with  the  common  notation 
for  stiffness  and  damping  coefficients,  may  be  expressed  as 

h  -  +  (2) 

where  is  the  stiffness  matrix  and  is  the  damping  matrix,  v  is  the 
frequency  of  vibration.  Most  commonly,  lateral  lineal  and  angular  displace¬ 
ments  do  not  interact  with  each  other  so  that  the  non-vanishing  portions  of 
and  are  separate  2x2  matrices.  That  is 


■N 


(£n> 


lineal 

0 

0 


0 

0 


(*N> 


0 
0 

angular  J 


(3) 


2.2  General  Bearing  Support  Characteristics 

The  function  of  a  bearing  is  to  restrict  the  rotor  axis  to  a  nominal  axis 
under  realistic  static  and  dynamic  load  environments.  Deviation  of  any 
particular  point  of  the  rotor  axis  from  the  nominal  line  can  be  character¬ 
ized  by  three  lineal  and  two  angular  displacements.  These  may  be  designated 

as(6,6,6,0,O)in  accordance  with  a  right-handed  Cartesian  reference 
x  y  z  x  y 

system.  The  z-coordinate  is  coincident  with  the  reference  axis  and  is 

directed  toward  the  spin  vector.  (0  ,  ©  )  are  rotor  axi9  inclinations 

x  y 

respectively  in  the  z-x  and  z-y  planes.  The  x-coordinate  is  directed  toward 
the  predominant  static  load;  e.g.,  earth  gravity.  Ideally,  the  bearing  would 
resist  the  occurrence  of  any  displacement  so  that  the  reaction  force  system 
imparted  by  the  bearing  to  the  rotor  is  generally  expressed  in  matrix  notation 
as 

F  *  Z  .  x  (10) 


F  is  a  column  vector  comprising  the  five  reaction  components  (F  ,  F  , 

x  y 


M  ,  M  ),  while  x  is  the  displacement  vector  (6,6,  6_,  0^,  0_)  . 
x  y  x  y 


a  (5  x  5)  matrix  containing  the  elements  Z 


ij 


,0,0).  Z  is 
z  x  y  ■ 

with  both  indices  (i,  j) 


ranging  from  1  to  5.  The  values  of  characterize  how  rotor  displacements 
are  being  resisted  by  the  bearing. 


From  the  standpoint  of  dynamic  perturbation,  distinction  is  made  between 
a  static  equilibrium  component  and  a  dynamic  perturbation  component  for 
both  the  displacements  and  the  reactions.  Thus, 

x  -  x  +  x’ ;  F  •  F  +  F'  (11) 

—  -o  —  —  -o  — 

(x' ,  £')  are  respectively  presumed  to  be  infinitesimal  in  comparison  with 
(xQ,  Fq) .  Accordingly,  Z^  are  regarded  as  dependent  on  xq  but  not  on  x' • 
To  illustrate  the  idea  of  perturbation  linearization,  one  may  examine  the 
one-dimensional  load-displacement  curve  shown  in  Figure  2. 


As  illustrated,  the  load-displacement  relationship  is  a  10/9  power  lav  in 
accordance  with  the  Hertzian  point  contact  formula.  It  is  not  possible 
to  describe  the  entire  range  by  a  linear  approximation.  However,  if  a 
small  dynamic  perturbation  is  taken  around  a  static  equilibrium  point, 

6'  <6  ,  the  small  segment  of  the  load-displacement  curve  can  be  approxi« 

X  Xq 

mated  by  a  local  tangent  line.  The  corresponding  force  increment  is 


F' 

x 


3F 

_ x  X. 

36  6x 
x 


(12) 


where  6'  is  the  incremental  displacement.  3F  / 3 6  will  depend  on  the 

X  XX 

amplitude  of  6 

o 

The  question  of  history  dependence  is  resolved  by  regarding  x'  as  periodic 
motions  at  any  frequency  v  of  interest,  and  Z accordingly  would  have 
both  real  and  imaginary  parts  and  may  also  be  dependent  on  both  the  rotor 
speed  u  and  the  vibration  frequency  v. 

To  avoid  notational  clumsiness,  the  primes  will  be  dropped  from  (F' ,  x') 
which  are  understood  to  be  dynamic  perturbation  quantitites  unless  the  sub¬ 
script  "o"  is  used  to  designate  the  static  equilibrium  condition. 


2.3  Cylindrical  Roller  Bearing  Characterization 

In  many  ways  the  cylindrical  roller  bearing  is  much  simpler  to  model  from  a 
rotor  dynamic  point  of  view  than  a  fluid  film  bearing.  In  general,  the 
following  two  simplifications  can  be  made: 

.  The  restraining  characteristics  do  not  include  cross  coupling 
among  the  various  displacemenc/inclination  coordinates. 

.  The  restreining  force/moment  is  normally  temporally  in  phase 
with  the  displacement/inclination. 
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Figure  3  illustrates  a  typical  radial  roller  bearing  haviug  cylindrical 
rollers.  In  the  bearing  illustrated  the  rollers  are  guided  by  closely 
fitting  flanges  on  the  outer  ring.  Other  constructions  are  common  and 
all  carry  load  in  the  manner  to  be  described. 


The  bearing  is  referred  to  an  orthogonal  xye  coordinate  system  with  the 
bearing  .  xis  directed  along  z.  The  bearing  reacts  to  the  linear  displace¬ 
ments  6  and  <5  and  to  the  rotations  or  raisalignmer  ts  G  and  0  .  All  are 
x  y  x  y 

shown  in  their  positive  senses  in  Figure  4. 


Figure  5  illustrates  the  roller  index  q. 


2.3.1  Stiffness 


The  total  characterization  of  a  cylindrical  roller  bearing's  stiffness 
can  be  expressed  by  the  matrix. 


[K] 


3F 

3F 

3F 

3F 

X 

X 

X 

X 

3x 

3y 

3z 

3G 

3F 
— Z 

3x 


3F 
— H 
3y 


3F 

9  7. 


3F 


3F 

3F 

3F 

z 

z 

z 

3x 

3y 

3z 

3F 

3F 

z 

z 

30 

30 

X 

y 

3M 

3H 

3M 

_ X 

_ X 

_ X 

3x 

3y 

3z 

3M 

X 

3M 

X 

30 

X 

30 

y 

3M 

3M 

3M 

3H 

_ )L 

V 

-JL 

3x 

3y 

3z 

30 

V 

(13) 
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Figure  3.  Cylindrical  Roller  Bearing 


Figure  4.  Bearini 


1 


The  lineal  and  angular  etlffneao  matrices  (Equations  5  and  7)  can 
be  derived  from  Equation  (13).  For  example: 


CK) 

lineal 


(K) 

angular 


3F 

3F 

X 

X 

3x 

3y 

3F 

3F 

— SL 

—L 

3x 

3y 

_3M 

3M 

X 

X 

30 

30 

X 

y 

3M 

3M 

—X 

30 

30 

X 

y 

(14) 


(15) 


Note  that  the  axial  components  of  stiffness  are  zero  for  a  cylindrical 
roller  bearing.  However,  these  terms  have  been  retained  on  the  general 
cylindrical  roller  bearing  stiffness  matrix  for  consistency  with  the 
ball  and  tapered  roller  bearing  analyses  included  in  earlier  volumes  of 
this  series  (4,6). 


2,3.2  Damping 

An  extensive  search  of  the  literature  revealed  no  experimental  damping 
data  for  cylindrical  roller  bearings.  As  the  current  state-of-the-art 
does  not  permit  accurate  calculation  of  cylindrical  roller  bearing 
damping,  no  damping  data  is  included  in  this  report. 


2.4  Cylindrical  Roller  Bearing  Under  Radial  Loading 

Solution  for  the  stiffness  matrix  of  a  cylindrical  roller  bearing  under 
radial  loading  is  a  tedious  problem  and  requires  the  use  of  a  digital  com¬ 
puter.  In  this  section,  the  derivation  of  the  solution  is  described.  A 
computer  program  for  obtaining  the  solution  is  included  in  the  Appendix. 
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2.4.1  Bearing  Applied  Forces  and  Moments 


In  Che  present  problem  an  external  force  Is  applied  to  the  inner  ring 
along  x,  only,  and  the  inner  ring  is  constrained  to  movement  along  x. 
There  may,  however,  be  initial  displacements  in  other  modes  which  are 
fixed  and  which  result  in  reactions  in  other  modes  as  will  be  shown. 

As  the  result  of  the  four  displacements  described  previously  in 

Figures  3  and  4,  there  are  the  reactions  F',  F',  M' ,  and  M'  .  F'  and 

x  y  x  y  x 

F'  are  forces.  M'  and  M'  are  moments.  All  are  shown  in  their  posi- 
y  x  y 

tive  sense  in  Figure  4.  External  force  F^  is  applied  at  the  inner 

ring  center.  The  senses  of  the  signs  are  the  same  aa  for  the  reaction 

F' . 
x 


2.4.2  Roller  Geometry 

Figure  6  shows  the  boundary  dimensions  of  a  typical  cylindrical  roller. 
Roller  mass,  moment  of  inertia,  and  location  of  the  center  of  gravity 
are  calculated  assuming  the  roller  is  a  flat-ended,  cylinder  bounded  by 
d  and 

I  ia  the  effective  length  of  the  roller  and  represents  the  maximum 
working  length  of  the  roller.  The  contact  pattern  between  roller  and 
race  must  lie  within  the  effective  length. 

is  the  total  length  of  the  roller. 

ip  is  the  length  of  the  flat  portion  of  the  roller  profile  and  may  be 
zero  for  a  fully-crowned  roller. 

Rc  is  the  roller  crown  radius  and  is  struck  from  the  central  plane  of 
the  roller.  V  is  the  drop  of  the  crown  measured  at  the  extremities  of 
the  effective  length.  Thus  we  can  write: 
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where 


R2  -  -£ 
c  2 


The  crown  drop,  V,  is  used  In  checking  the  roller  profile.  If  the 
crown  radius  is  known,  the  drop  at  a  distance  G  from  the  extremity 
of  the  effective  length  is 

R2.(ii.0]2 


If  the  drop  V'  at  G  is  known  the  crown  radius  is 


*  2  £  2  ; 

— ~  »  G  —  —  -  p' 

2  12 


B  is  the  corner  break  which  is  the  same  at  both  ends  of  the  roll.  The 
exact  form  of  the  corner  break  is  unimportant  as  long  as  the  corner 
blends  smoothly  into  the  roller  crown.  The  effects  of  the  corner 
breaks  are  neglected  in  calculating  roller  mass.  The  effect  of  the 
crown  is  also  neglected. 

The  mas9  of  the  roller  is 
U  +2B)d2p 

m  “  491.98  (20) 

3 

where  p  is  the  material  density  of  the  roller  in  lb/ in  . 

The  centrifugal  force  acting  on  the  roller  is 


E  is  the  pitch  diameter  of  the  bearing  in  inches, 
velocity  of  the  roller  and  is 


ft  is  the  orbital 


ni(l+T)  +  n2  (1-Y) 


where 


Y  - 


d 

E 


(23) 


and  n£  are  Che  angular  velocities  of  outer  and  inner  rings,  respectively, 
in  radians/sec. 

The  subscripts  1  and  2,  used  in  Equation  7  and  subsequently,  refer  to  outer 
and  inner  contacts  respectively. 

Figure  7  shows  the  forces  and  moments  acting  on  a  roller  which  is  in 
contact  with  both  outer  and  inner  races. 

and  P2  are  the  contact  loads.  and  M,,  are  the  contact  moments 
resulting  from  non-uniform  loading  along  the  roller's  length.  F£  is  the 
centrifugal  force  acting  at  the  mass  center  of  the  roller. 

In  the  present  problem  the  radial  roller  bearing  is  loaded  by  a  single 
force  applied  along  x  and  working  displacements  are  limited  to  a  single 
deflection  along  x. 

Assume,  for  the  moment,  that  working  displacements  are  possible  along  x 
and  y  and  about  x  and  y.  Then  the  approach  of  inner  and  outer  race  for 
a  roller  at  azimuth  *  is 

P 

A  “  (6  +6")  cos*  +  (6  +5")  sin*  -  ■=—  (24) 

xx  y  y  2 

where  PD  is  the  diametral  clearance. 

The  relative  misalignment  of  the  inner  race  with  respect  to  the  outer  is 

0  -  (0  +0")  sin*  +  (0+0")  cos*  (25) 

xx  y  y 
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4 


i 


Forces  and  ?toments  Acting  on  Roller 
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The  double  primed  Items  in  Equations  24  and  25  are  initial  displacements. 
If  initial  displacements  exist  for  a  particular  mode,  working  displace- 
ments  in  that  mode  are  prevented. 


For  the  present  problem  the  linear  displacement  6^  and  the  rotations  0^ 
and  0y  are  non-existent  and  the  approach  and  misalignment  of  the  races 
are,  for  this  case, 


4  -  (6x+6Jp  cos$  +  Sy  9in*  "  T" 


(26) 


and 


0-0^  sin*  +  0^  cos*  (27) 

When  0  is  positive  the  left  end  of  the  roller,  viewed  with  the  axis 
horizontal,  tends  to  be  pinched. 


If  A^  is  the  approach  of  the  roller  to  the  outer  race;  the  approach  of  the 
roller  to  the  inner  race,  at  the  midpoint  of  the  roller  is 

A2  -  A-Ax  (28) 

and  the  misalignment  at  the  inner  contact  is 

©2  -  O-Oj^  (29) 

Figure  8  illustrates  the  geometric  intersection  of  a  roller  and  raceway. 

The  profiles  of  race  and  roller  bodies  are  referred  to  an  XY  coordinate 
systen.  Note  that  the  X  axis  is  positive  to  the  left  of  the  origin. 

The  equation  of  the  race  surface  is 


Y  - 


0 


(30) 


The  equation  of  the  flat  portion  of  the  roller  is 


Y  -  A  +  XtanOi 


(31) 


The  equation  of  the  crowned  portion  of  the  roller  is 


(X  -  Hsin©.)2  +  (Y  +  H  cos© .  -  A,)2  -  R2 

l  lie 


(32) 


In  the  above  the  subscript  i  is  1  for  an  outer  contact  and  2  for  an 
inner  contact. 


The  intersections  of  the  race  and  the  crowned  surface  of  the  roller  occur 


at  X.  and  Xg 


i  m~\jK  - 


+  HainS. 


(33) 


R2  -  (HcosQ .-A.) 2  +  HsinO. 
c  li  l 


(34) 


X  and  x  must  be  within  the  projected  extremities  of  the  roller 
Ai  Bi 


crown. 


That  is 


X.  $  X. 
Ai  Ai 


\ 
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1 


where 


.«*  ..  .  ' 


★ 

X,  *»  -r~  cos0.  +  Vsin©. 
A.  2  l  i 


B . 


(37) 


—  cos0^  +  7sin0^ 


(38) 


If  the  quantity  under  the  radical  in  Equations  33  and  34  is  zero  or 
negative  there  is  no  contact  between  roll  and  race. 


If  tt~  cosO.  £  X.  there  is  also  no  contact. 
L  i  A^ 


If  X.  >  X  ,  X  is  set  equal  to  X.  . 
i  Ai  Ai  Ai 


If  Xg  <  Xg  ,  Xg  is  set  equal  to  Xg  . 


If  —  cosQ^  >  Xg  >  -  —  cos©i  and 

X,  >  -r—  cosO.  the  value  of  X„  i9 
Af  2  i  ^B1 


Y 


tan©. 


39) 


From  Figure  7  the  conditions  for  roller  equilibrium  are 


-  P,  +  P,  +  F  -  0 
12c 


(40) 


M  -  M2  -  0 


(41) 


Equations  40  and  41  are  a  set  of  non-linear  simultaneous  equations  in 
which  the  variables  are  A^  end  0^  at  the  outer  contact  of  the  particular 
roller  under  consideration. 


i 

i 
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..  ^  ?  , r> ^wn?*** yf  *i*m**r+> 


r  ^ 


•-_  «s..  ^  '  -f  \ 

•^jfc'  *#.  L  ""'  T. 


From  Figure  8,  the  intrusion  of  the  roller  into  the  raceway  is 


A  ■  A.  +  Xtan0. 
x  i  i 


|X|  $  y  cos©t 


f  2  o 

R*  -  (X-Hsin©)  -  Hcos©i  +  |x|  >  ~  cos©1  (43) 


The  derivatives  of  with  respect  to  ©^  will  be  required  later. 
They  are 


d©.  2_ 

i  cos  0, 


|X|  5  f-  cos ©1 


(X-Hsir.Q^)  Hcos©^ 

R2  -  (X-Hsin© . )2 
c  i 


+  HsinQ . 


|x|  >  f  cos©i 


Lundberg  (7)  gives  the  approach  Ax  of  two  cylindrical  bodies  pressed 
together  with  the  uniform  loading  px  as 

(n  +  nE)  f  |  X.-X_r 

Ax  -  — 5T“  Px  1-8864  +  in  (46) 

'  x 

_  J 

nR  and  nE  are  elastic,  constants  for  race  and  roller,  or  respectively,  having 
the  form 

(l-vR2)  <l-vE2) 

nR-— °*nE-  — —  (47) 

where  v  is  Poisson's  ratio  and  E0  „  is  the  modulus  of  elasticity. 

R,  L 

b  is  the  semi-width  of  the  pressure  area  in  the  rolling  direction 


7.  G.  Lundberg:  Elastische  Beruhrung  ziveier  Halbraume,  VDI  Forschung, 
Sept. /Oct.,  1939. 
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(48) 


V 


P  d  (1  +  C 


iYi>] 


1/2 


is  1  for  an  outer  contact  and  -1  for  an  inner  contact.  Y  is  given 
by  Equation  23. 


The  value  of  p  corresponding  to  A  is  required.  It  cannot  be  obtained 

X  X 

from  Equation  46  in  closed  form.  It  can  be  obtained  numerically  as 
follows. 


Let  p'  be  an  estimate  of  p  .  A  good  starting  value  is 

*  moV>«  * 

.  x 


(vv 


1/9 


(49) 


An  improved  value  of  p^  is 


Px  * 


(A'-A  ) 

X  X 

dP 

_ x 

dp' 

X 


(50) 


A^  is  the  approach  of  race  and  roller  bodies  calculated  for  the  current 
estimate  p^  using  Equation  46. 

dA' 

-t-t  is  obtained  from  Equations  46  and  48  using  the  current  estimate  p' 
anS  is 


dA' 

K 


(W 

2tt 


1.3864  +  In 


Vxb1 


2b 


(51) 


Iteration  of  Equation  50  yields  p^  to  any  desired  accuracy. 
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f 


The  contact 


£orce  P  and  the  moment  M  are 


P 


1 


M 


i 


p  dX 
*x 


p  dX 
rx 


(52) 


(53) 


Equations  40  and  41  may  now  be  solved  numerically  for  A^  and  0^. 


If  estimates  are  made  of  the  variables  A^  and  P^  Equations  40  and  41 
may  not  be  satisfied  and  there  will  be  the  residues  and  £2  for 
Equations  40  and  41  respectively.  Differentiating  Equations  40  and 
41  gives: 


del 

dPl 

4. 

dP2 

dA2 

dA,  ' 

~dAl 

dA2 

dAl 

de^ 

dPl 

_i_ 

dP2 

do2 

dP^  “ 

-  dGx 

T 

d°2 

d°l 

dc2 

dMx 

dM£ 

dA2 

dAl  " 

dA2 

dAl 

de2 

dHL 

dM2 

d°2 

d0l 

d0^ 

d0l 

where 


(54) 

(55) 

(56) 

(57) 

(58) 
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If  A^  and  0|  are  current  estimates  of  the  variables,  improved  values 
are  A^  and  0^ 


AI- 


3i- 


de. 


d^ 

de. 


d0. 


de, 


(60) 


d*l 

d3l 

de2 

de2 

dAi 

d  Q1 

d£l 

dA, 

E1 

de2 

dAx 

e2 

111 

del 

dA^ 

d©x 

de2 

d£l 

dAl 

d0l 

(61) 


The  right  members  of  Equations  60  and  61  are  evaluated  at  current 
estimates. 


Iteration  of  Equations  60  and  61  yields  A,  and  r. 


to  any  desired  accuracy 


The  derivative*  of  and  with  respect  to  A^  and  are 


dpi 

aa 

Ai 

d\ 

d*x 

^  -  : 

k 

dA 

X 

dAi 

dpi  \ 

'XA 

Ai 

i 

dp 

rx 

dA 

X 

dei  -  < 

k 

dA 

X 

d0i 

(62) 


(63) 


If  Equations  27,  28,  40  and  41  are  differentiated  with  respect  to  A, 
there  results  four  simultaneous  equations  which  are  linear  in 

ddL  d A2  dOx  d©2 

•77—,  -rr~t  TT~»  and  7—  and  from  which  the  latter  can  be  obtained. 
dA  dA  dA  dA 

dAL  dOj 

Only  and  are  of  interest  here. 
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i 

*• 


% 


dpl 

ib.  + 

dP2 

dAl 

dA  + 

dA2 

dMx 

dAi 

dM2 

d^ 

dA 

dA2 

dAi 

+  ^2 

dA 

dA 

fl 

d^ 

a  4 

dA 

dA 

*  ( 

Equations  66  through  69 


dA, 

dP, 

de. 

dP 

2 

1 

i 

dA 

d<3, 

dA 

+  de! 

dA, 

dM, 

de. 

dM 

2 

+  — i 

l 

dA 

d9l 

dA 

dO 

d  A 

easily  solved  for  — 


]<«!  f-; 


t  *:  ■r'-. ' - - - - 


The  conditions  for  equilibrium  under  the  single  radial  load  F  is 


F'  +  F  -  0 
x  x 


The  only  working  response  of  the  bearing  is  the  linear  displacement  6^ 


If  6'  is  an  estimate  of  6  an  improved  value  of  6  is 

X  X  X 

<f;  +  f  > 

»  #  f  A  A 


The  right  member  of  Equation  76  is  calculated  at  current  values  of  6 

Iteration  of  Equation  76  yields  6*  to  any  desired  accuracy. 
dF' 

Although  only  j~-  is  required  in  the  foregoing,  all  the  derivatives 
x 

are  required  to  complete  the  matrix  needed  for  stiffness  calculations. 
The  elements  of  the  complete  matrix  are 


dFx  5T*  q 

,5  ,0  ,o  )  *  c°8*q  d<{v*6v*°v»0v) 

x  y  x  y  x  y  x  y 


x’  y'  x  y 


u 

Y'  sin$ 


d(6x’V0x’V  H  d(6x,6y,0x,0y) 

q*l 


(80) 


dM' 

x 


d(5x’VVV 


n  dMl 

Z.  “’'iJvTv5? 


dNy  V'  dMlq 

4<4.>VVV  *  L  co,*q  d(5x,«y,ex.ey) 


(81) 


d(F' ,F' ,F ' ,M' ,M' 
x  y  z  x  y 

d6 


(82) 


where 


dP, 


d(6x,«y,0x.°y) 


r  dPj_  dAj_ 


d A.  d A 

K  d 


dPj_  d01  -i 


d0,  d  A 

1  q 

n  n 


dA 


d(<$  ,6  ,0  ,0  ) 
x  y  x  y 


\ 

i  dA,  dO 

L  XQ  q 


d?1  d01  1 


dO,  d0 
1  d 


d0 


d(6  ,6  ,0  ,0  ) 
x  y  x  y 


(83) 


dM, 


dM,  dA, 


r 


d(6x,ay,6x,ey)  L  dAlq  dAc 


dMx  d0^  "i 


d0,  dA 

K  d 


dA 


d(6  ,6  ,0  ,0  ) 
x  y  x  y 


r  dMx  dAx 


dA,  d0 
1  d 


dM^  d0^  i 


d0,  d0 
1-  d 


dO, 


d(5  ,5  ,0  ,0  ) 
x  y  x  y 


(84) 


The  derivatives  of  A^  and  0^  with  respect  to  the  inner-ring  displacements 
are,  from  Equations  24  and  25. 
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d(e  ,0  ) 

x  y 


% _ 

d«x’V 


(88) 


d0 


d0 


(89) 

(90) 


Some  rollers  may  be  out-of-contact  with  the  inner  race  and  are  loaded 
against  the  outer  race  by  centrifugal  force.  Such  rollers  contribute 
to  the  bearing's  reactions  but  have  no  effect  on  the  derivatives  of 
P  and  M.  . 

q  q 

The  complete  matrix  for  stiffness  calculations  is  then: 


As  discussed  earlier,  the  axial  components  of  stiffness  are  zero  for 
a  cylindrical  roller  bearing. 
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SECTION  III 

APPLICATION  OF  COMPUTER  PROGRAM 

The  analysis  of  Section  II  has  been  programmed  in  Fortran  IV  for  a 
digital  computer  and  is  suitable  for  use  on  the  CDC  6600.  A  program 
listing  is  presented  in  the  Appendix. 

3.1  Sample  Test  Case 

To  illustrate  a  typical  case  consider  the  cylindrical  roller  bearing  in 
Figure  9.  This  is  a  cylindrical  bearing  typical  of  a  design  employed 
in  high  speed  applications.  The  geometry  of  this  sample  bearing  is 
summarized  below. 


Number  of  rollers,  n 

Roller  diameter,  d 

Pitch  diameter,  E 

Total  length  of  roller ,  £T 

Length  of  flat  portion  of  roller  t F 

Roller  crown  radius 

Roller  corner  break 

Crown  drop  gage  point 


.6299  in. 
7.75  in. 
.6299  in. 
.2  in. 

17  in. 

.05  in 


The  operating  conditions  for  the  sample  case  are: 

Diametrical  clearance,  PD  *  .0058  in. 

Rotational  speed,  N^  ■  13,230  rpm 

Radial  load,  F  ■  1,465  lbs. 
x 

3.2  Input  Format 

Figure  10  presents  the  input  data  format  and  Figure  11  shows  the  actual 
input  data  for  the  sample  case  operating  conditions. 


3.3  Output  Format 

Figure  12  presents  the  output  data  for  the  sample  case.  The  input  data 
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Pitch  Diameter 
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Problem  Data  Input 
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Figure  12  Output  Data  for  Sample  Problem 
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Figure  12  Output  Data  for  Sample  Problem  (Continued) 


are  summarized  In  Figure  12,  followed  by  the  output  data  Including  the 
internal  load  distribution  aa  well  as  various  other  stress  and  displace¬ 
ment  parameters.  The  stiffness  matrix  is  given  on  the  last  page  of 
Figure  12. 
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APPENDIX 


COIffUTETl  PROGRAM 
FOR 

CALCULATING  STIFFNESS  MATRIX 
OF 

CYLINDRICAL  ROLLER  BEARING 


i 

I 

\ 
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